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—.^ I Abstract. We discuss the implementation and validation of a procedure 

Q ■ aimed to detect groups and clusters of galaxies in DPOSS catalogs. 

o' 
o 

Qh' !• Introduction 

;_i ■ The Digitised Palomar Sky Survey (hereafter DPOSS) offers an unique data set 

c/3 ! to explore the large scale structure of the nearby {z < 0.2) universe due to its 

^ ■ virtually unlimited sky coverage, to its deepness (21.5 in g) and the availability 

of three colours (J, F and N). 

In this poster we discuss a first application, on a small part of the whole data set, 

of a newly implemented procedure to search for candidate galaxy clusters and 
C^ ' groups. Final goal of this project is a robust measurement of the multiplicity 

function (hereafter MF) of galaxies. 

The MF is a powerful tool to test the various cosmological scenarios. For istance, 

if the protostructures form from a hierarchical gravitational instability caused 

by primordial density fluctuations: 

^ocA- m(-i/2+"/6) 

Pm 

where n is the initial spectrum index at era of baryonic recombination and if one 
identifies groups of galaxies with mean density rg > k times the average cosmo- 
logical density pM, then the function of the mass distribution doesn't depend 
on the cosmic density parameter Q but only on the spectrum index n (Gott & 
Turner 1977). The groups which are seen now at a density enhancement k would 
therefore correspond to density fluctuations of amplitude A at recombination. 
Therefore, a measurement of the MF at the present epoch derived from a cata- 
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logue of groups of galaxies identified for a given k, gives a measure of the initial 
spectrum index n. 



2. Detection of candidate groups 

Even though on rather arbitrary ground, we define as a group all galaxy ensem- 
bles having less then 15 members having total magnitude within 3 magnitudes 
from that of the brightest members. The implemented procedure is described 
in detail in De Fihppis (1999). 

In order to compile a catalogue of candidate groups of galaxies in absence of 
redshift information, we implemented a modified version of the van Albada al- 
gorithm (Soares D.S.L. 1989). This algorithm makes use of apparent magnitudes 
and projected positions on the sky only, and gives for each pair of adjacent galax- 
ies their probability of being physically bound. 

Assuming a Poisson statistic, the probability that the distance of a fixed galaxy 
to its nearest non-physical companion lies between and Q + dQ, is given by: 

Pi{e)de = exp[-7:e'^n]27:ende 

The introduction of an adimensional distance x allows to combine the angu- 
lar separation of different pairs to a single distribution removing the effects of 
clustering in the background. 



3. Detection of candidate clusters 

Starting points are the individual J, F and N catalogues produced by SKICAT 
for a given DPOSS field after conversion to the g, r and i Gunn-Thuan system. 
The three catalogues are then matched (assuming a maximum matching distance 
of 7", for details see Puddu et al. 1999). 

After matching, all objects having r > 19.5 are filtered out and the resulting 
matched catalogue is binned into equal area square bins of 1.2'to produce a 
density map. Then, S-Extractor (Bertin & Arnouts 1996) is run in order to 
identify and extract all the overdensities having number density 2a above the 
mean background and covering a minimum detecting area of 4 pixs (4.8'sq.). In 
this way, as in the Schectman (1985) approach, we are not assuming any a priori 
cluster model. An application to the DPOSS field 610 is shown in Figure 1. 
All the previously known Abell and Zwicky clusters are recovered and many 
new candidates detected. For each candidate cluster we then measure the S/N 
detection ratio and the Abell richness class. 



4. Validation of the algorithms 

In absence of a suitably complete samples of galaxies with known redshift we 
were forced to validate the algorithms shortly summarised above using either 
simulations or photometric techniques. In what follows we shall discuss first the 
groups and then the clusters. 
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Figure 1. Results for DPOSS field n. 610 (1^+15°). Circles mark 
the position of detected candidate clusters while squares mark already 
known Abell and Zwicky clusters. 

4.1. Group algorithm validation 

In order to test the accuracy (lost groups) and reliability (spurious groups) of 
the algorithm we tested it on 70 realistically simulated fields. We first produced 
the galaxy background by assuming uniform distribution and the field luminos- 
ity function given by Metcalfe et al. 1995; then we added groups of galaxies 
according to the multiplicity function by Turner E.L. &: Gott J.R. 1976, with 
redshift computed according to: 
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and absolute magnitude of the brightest galaxy in the group taken from the 
cumulative LF for groups of galaxies: 

^M)dM = $*[ioO-4(A^*-A^)]("+i) . exp[-10°-4(*^*-*^)]dM 

where M* = —20.85 and a = —0.83 — 1-0.17. Other parameters varied in the 
course of the simulations were: i) the diameter of the group inside a gaussian 
distribution centered at Dq = 0.26Mpc; ii) the maximum possible redshift for 
a group (chosen in the range 0.2 < z < 0.7). Different simulations have been 
performed changing the lower limit of the probability for wich two galaxies are 
considered physical companions by the algorithm. The best compromise be- 
tween accuracy and reliability is found at p = 0.6. 
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Figure 2. Comparison of our MF against those derived by other au- 
thors. 



In Figure 2 we show the multiphcity function obtained from 12 DPOSS plates 
(covering a total sohd angle of ASSdeg"^ ) compared to six other MFs taken from 
the hterature: MF derived from magnitude hmited surveys making use of red- 
shift information, the CfA survey (Geher & Huchra 1983; Garcia et al. 1993; 
Ramella &: Pisani 1997; Ramella et al. 1998), MF derived from diameter limited 
survey (Maia & Da Costa 1989) and MF derived from limited magnitude survey 
without redshift informations (Turner & Gott 1976). 

As it can be clearly seen, the MF derived by appling our method to DPOSS ma- 
terial is at 98% confidence level undistinguishable from literature MF obtained 
from redshift surveys. 



4.2. Cluster algorithm validation 

We took two different approaches: through simulation and through magnitude 
- color diagrams. In order to produce realistic simulations we followed a rather 
lengthy procedure. First of all, we needed a fairly accurate multiplicity function 
to be used as input parameter for the simulation. Since there is no good MF 
available in literature, we used a set of 20 DPOSS triplets to identify clusters 
and derive a preliminary MF (we used only very high S/N clusters for which 
there is no doubt about their physical nature). 

In order to simulate the fields we first distribute randomly the cluster centers 
and we assume the cluster distribution to be uniform in the volume element and 
randomly assign to each cluster a redshift in the range [0.02, 0.2]. Then we grow 
the clusters around their centers randomly selecting a radius in the range [2 — 3] 
Mpc and weighting the number of members accordingly to the above derived 
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MF. Absolute magnitudes for galaxies are then extracted from a Schechter func- 
tion (M* = —21.41, a = —1.24) and scaled to apparent magnitudes taking into 
account the K-correction term. We finally add field galaxies using the Metcalfe 
field counts and LF. Even though simulations are still in progress, preliminary 
results show that our algorithm can recover 95% of the clusters at z < 0.15 and 
70% of the clusters in the range [0.15,0.2]. 



4.3. Photometric validation 

As stressed above, DPOSS data are calibrated via independent CCD frames 
taken at various telescopes in the 1 — 2m range. This means that we have at our 
disposal a large number of deep fields, where galaxy photometric properties can 
be measured with high photometric accuracy. In the last two years we selected 
calibration fields in order to largely overlap with our candidate clusters sample. 
These data were used to derive color - magnitude diagrams. 
To calibrate the CCD frames taken in the g, r and i Gunn-Thuan system, we 
use S-Extractor to identify galaxies and derive the colour-magnitude diagram (r 
vs. g — r) for the candidate cluster region and for a test field (background) well 
selected. Then we statistically subtract the field contamination and identify the 
excess of galaxies as cluster members. The existence of an excess of galaxies in 
the putative cluster line of sight and the fact that they obey to the usual colour- 
magnitude relation is taken as confirmation of the existence of the putative 
cluster. These galaxies are then used to build the cluster radial profile. The 
colour-magnitude diagram allows also to derive a rough photometric redshift 
(which, due to the broad band, has however a rather large error). 
Preliminary results confirm that for z < 0.15 the candidate clusters found by 
our algorithm are 95% "true" clusters. At higher redshift (0.15 < z < 0.2) this 
fraction drops to 70%. 

We plan to improve the simulation by adding non-circular simmetry for 
the simulated cluster and to study the dependence on the zero-order correla- 
tion functions. We also plan to apply the colour- magnitude diagram test to all 
clusters avalaible in our calibration data set. 



References 

Bertin, E., Arnouts, S. 1996, A&AS, 117, 393 

De Filippis, E., Lo ngo, C, Andreon, S., et al. 1999, MmSAI, in press ( |astro- 
ph/990936q ) 



Garcia, A.M. 1993, A&AS, 100, 47 

Geller, M.J., Huchra, J.P. 1983, ApJS, 52, 61 

Gott J.R., Turner E.L. 1977, ApJ216, 357 

Maia, M.A.G., Da Costa, L.N. 1989, ApJS, 69, 809 

Maia, M.A.G., Da Costa, L.N. 1990, ApJ, 352, 457 

Metcalfe, N., Shanks, T., Fong, R., Roche,N. 1995, MNRAS, 273, 257 

Puddu, E., et al. 1999, MmSAI, in press ( ^stro-ph/9909367] ) 

Schectman, S.A. 1985, ApJS, 57, 77 



E. Puddu et al. 









1 


• *Ha 


■•■\-^ : 




■■■'-%\ 








.1 field!.. ^ -- 








20 2a 


• • > . 


Jb. 






,__,,,,,,,,.,_ 



20 22 



, .'■'^W>■•■ 




Figure 3. Color - magnitude diagram for the candidate cluster ob- 
tained accordingly to the procedure described in the text. Right (high 
S/N candidate): the cluster turns out to be of richness class I and to 
be at z = 0.10 — 0.12. Left (low S/N candidate): the cluster turns out 
to be of richness class I and to be at z = 0.26 — 0.3. 



Soares, D.S.L. 1989, PhD Thesis 

Turner E.L., Gott J.R. 1976, ApJS, 32, 409 

Weir, N., Fayyad, U.M., Djorgovski, S.G., Roden J. 1995, PASP, 107, 1243 



